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In the quantum anomalous Hall effect, chiral edge modes are expected to conduct spin polarized 
current without dissipation and thus hold great promise for future electronics and spintronics with 
low energy consumption. However, spin polarization of chiral edge modes has never been established 
in experiments. In this work, we theoretically study spin polarization of chiral edge modes in the 
quantum anomalous Hall effect, based on both the effective model and more realistic tight-binding 
model constructed from the first principles calculations. We find that spin polarization can be 
manipulated by tuning either a local gate voltage or the Fermi energy. We also propose to extract 
spin information of chiral edge modes by contacting the quantum anomalous Hall insulator to a 
ferromagnetic (FM) lead. The establishment of spin polarization of chiral edge modes, as well as 
the manipulation and detection in a fully electrical manner, will pave the way to the applications 
of the quantum anomalous Hall effect in spintronics. 

PACS numbers: 73.20.-r, 72.25.-b, 85.75.-d 


Introduction - For a two dimensional electron gas un¬ 
der a strong magnetic field, Landau levels are formed 
and drive the system into a state characterized by the 

zero longitudinal resistance and the quantized Hall con- 
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ductance with an integer value of This phenomenon 
is known as the quantum Hall (Qffi effect, which was dis¬ 
covered by Von. KlitziM in 1980 [l|. Recently, it was the¬ 
oretically predicted that this type of quantization in 
the Hall conductance can also be realized in magnetic in¬ 
sulating materials at a zero external magnetic field. This 
phenomenon, dubbed as the “quantum anomalous Hall 
(QAH) effect”, was soon observed experimentally in mag¬ 
netically doped topological insulators (TIs), the Cr or V 
doped (Bi,Sb) 2 Te 3 films [1-0. The physical origin of 
the QAH effect is spin-orbit coupling and exchange cou¬ 
pling between magnetic moments and electron spins in 
magnetic materials, rather than magnetic fields and the 
associated Landau levels Q. The experimental observa¬ 
tion of the exact quantization of Hall conductance and 
neglegible longitudinal resistance confirm the dissipation¬ 
less nature of transport for the QAH effect [H, 0 and the 
mapping of global phase diagram establishes the topo¬ 
logical eoniyalence between the QH effect and the QAH 
effect [^. flQ|. 

Similar to the case of the QH effect, dissipationless 
currents in the QAH insulators are carried by one dimen¬ 
sional (ID) chiral edge modes (OEMs), which propagate 
along one direction at the edge of the system and are 
robust against disorder scatterings. OEMs are believed 
to hold great promise for the potential applications in 
electronics and spintronics with low power consumption 
ini- Eor any spintronic application, it is required for 
OEMs to carry spin polarization (SP). Naively, one may 
expect that SP of OEMs should exist and follows bulk 
magnetization, but this issue has seldom been studied 
theoretically. In Ref. [l2|, SP of OEMs was studied in 


the context of a two band model, which is more relevant 
to cold atom systems. Eor condensed matter systems, 
it is more complicated since spin and orbital are mixed 
with each other due to spin-orbit coupling 

In this letter, we investigate SP of OEMs of the QAH 
effect in magnetically doped (Bi,Sb) 2 Te 3 . Surprisingly, 
we find that SP of OEMs though exists but does not fol¬ 
low bulk magnetization, and sensitively depends on the 
boundary conditions. In particular, we find that the di¬ 
rection of SP of OEMs can be manipulated by a local 
gate voltage, thus opening the possibility of controlling 
SP of OEMs in a fully electrical manner. We provide a 
simple physical picture of SP of OEMs based on the ef¬ 
fective four-band model and further study its behavior in 
the more realistic calculations based on the tight-binding 
model constructed by the Wannier function method of 
the first principles calculations 13 [m. We propose the 
spin valve effect [l^ of OEMs in a standard experimental 
setup of the Hall measurement with ferromagnetic (EM) 
leads to extract spin information of OEMs. 

Spin polarization of chiral edge modes - To study SP 
of OEMs in a QAH insulator, we first consider an effec¬ 
tive four band model for magnetically doped (Bi,Sb) 2 Te 3 
films 0 . The low energy physics of this system is domi¬ 
nated by two surface states from top and bottom surfaces, 
which hybridize with each other to open a gap due to the 
finite size effect. Meanwhile, the exchange coupling of 
electron spin arises from the doping of magnetic atoms. 
Thus, the effective Hamiltonian i^e// is given by 

Heff = VTz 0 {ky(Jx - kxCTy) + m{k)rx 0 (Tq pM ■ tq 0 a 

PV{x, y)Tz (8) ctq + A{k)ro (g) ao 

Ph{k^ P k^)To (g) (Tz, (I) 

where (Tx,y,z and Tx^y^z matrices are Pauli matrices of spin 
and layer (top or bottom) degree of freedom, and cfq = 
^0 = ^ 2 x 2 - m(k) = mo Pm 2 {k‘^ + ky) gives hybridization 
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between the top and bottom surface states. The M term 
describes exchange coupling between electron spin and 
magnetization of magnetic doping. We only consider the 
out-of-plane magnetization . A{k) = AQ + Ai{k‘^-\-ky) 
is a higher order correction of the surface states and h is 
defined as the strength of hexagon warping effect [illii. 
The V term denotes the asymmetric potential between 
the top and bottom layers which has spatial dependence 
and can be induced by a global or local gate voltage. 

The edge spectrum of our effective model, as well as the 
corresponding SP, can be evaluated with the help of the 
iterative Green’s function method [3113. We consider 
a semi-infinite system with the x direction still trans- 
lationally invariant and one edge parallel to x direction 
(x-edge). The following set of parameters, v = = 

0.5,mo = 0.1,m2 = 0.25, Aq = 0.15, = 0.05,h = 0.1, 

is chosen to keep the system in the QAH regime. The lo¬ 
cal Green’s function G{kx^uj) on the x-edge can be eval¬ 
uated iteratively. The total local density of states (DOS) 
and the spin DOS along the ith direction {i = x, y, z) are 
defined as po = — ;^Im[tr(G)] and pcr^ = —^lm[tr{G(7i)\^ 
respectively. As shown in the Supplementary Materials 
[l3|, the GEM can be easily identified from the peak of 
local DOS Po on the x-edge. The corresponding SP Si 
along the ith direction {i = x^y^z) for the GEM is nor¬ 
malized by the total DOS as Si = pajpo- In Eig. [T] 
(a) and (b), the SP of the OEMs is plotted for all three 
spin components with different local gate voltages and 
chemical potentials. Let us take the lattice constant to 
be unity and label the lattice site with an integer index 
n > 1, with (n — 1) being the distance between the nth 
lattice site and the boundary. Here the local gate voltage 
V is added only on all n = 1 lattice sites, which leaves 
the bulk states unchanged. We find that for a zero gate 
voltage (V = 0), only the 2 ; component (Sz) SP of OEMs 
is non-zero. If we apply a local gate voltage, SP becomes 
non-zero for both the y and z direction, but still keeps 
zero for the x direction. Therefore, SP of OEMs can exist 
within the y-z plane and can be controlled by a local gate 
voltage. It is also interesting to notice that the local gate 
voltage mainly controls the amplitude of Sy^ but barely 
change that of Sz (see Eig. [irb))[il- The chemical po¬ 
tential can also tune the magnitude of Sz (both Sy and 
5';^) at a zero (finite) local gate voltage, as shown in Eig. 
[1] (a). Therefore, our effective model (Eq. [T]) suggests 
that SP of GEMs can be manipulated by either applying 
a local gate voltage or tuning the chemical potential. 

Next we provide an analytical solution of the eigen 
wave function for the Hamiltonian (Eq. [T]) with V = 
Aq = Ai = h = 0 to understand the electrical tunability 
of SP of OEMs. Assume mo ,2 > 0, the system will enter 
QAHE regime when \Mz\ > mo- The zero mode of the 
Hamiltonian i^e// localizing on the edge can be solved 
exactly 0 as 

^{y) = C{e-^-+y-e-^-y)[\t)®{\ t,)) + |6)®(| I,))], (2) 




0.25 

EpfeV) 


- 0.2 


0 0.2 
VgateC^V) 


FIG. 1. Spin polarization Sx,y,z of chiral edge state in our 
effective model is plotted at fixed local gate voltage {V = 0.1 
eV) in (a), and at fixed Fermi energy {Ef = 0.35 eV) in (b). 




where G is a normalization factor and A± 

^/v‘^ 4m2(mo — Mz)). Here | t?/ (iy)) denote spin up 
(down) state along the y direction, and \t{b)) denotes the 
contribution from the top (bottom) layer of thin films. 
The wave function (j2j) of GEMs consists of two parts: one 
part on the top surface with SP along the -^y direction 
and the other on the bottom surface with SP along the 
—y direction. Thus, the SP is locked to the layer (top 
or bottom) for the OEMs. A local gate voltage can push 
the wave functions of GEMs into one layer, thus leading 
to the change of SP. 

The above analysis of SP of GEMs is based on the 
effective four-band model (Eq. [T]) and one may ask if 
these results still hold for a realistic system, such as Gr 
or V doped (Bi,Sb) 2 Te 3 . To answer this question, we 
carry out explicit calculations for a magnetically doped 
Sb 2 Te 3 thin film system with the realistic tight-binding 
model constructed from the maximal localized Wannier 
function (MLWE) method 14, TBI based on the first prin¬ 
ciples calculations, which has been successfully applied 
to the quantitative study of the QAH effect in magneti¬ 
cally doped (Bi,Sb) 2 Te 3 [2l|, [22[. The exchange coupling 
between electron spin and magnetization is included in 
the tight-binding model in the mean field approximation 
Hex = Here we consider only the con¬ 
tribution from out-of-plane magnetization and is 

the z directional spin operator for Sb (Te) atoms. In 
our calculation, we consider a film with two quintuple 
layers and choose the exchange coupling strength to be 
^Sh = 0-4 eV and Xtc = 0-0 eV, which is strong enough 
to drive the system into the QAH phase. The edge dis¬ 
persion is also calculated with the iterative Green func¬ 
tion method in a semi-infinite configuration along (100) 
direction, as shown in Eig. [2] (a). Eor edge modes, we 
find three left movers (the modes I, H and IV) and two 
right movers (the modes HI and V), suggesting that one 
chiral edge mode (left mover) and the other two trivial 
ID edge modes (or non-chiral edge modes) [2l[. Since 
the mode IV is directly connected to V, these two modes 
must be trivial ID edge modes and thus we focus on the 
modes I, H and HI below. The influence of the local gate 
voltage {V = 0.1 eV) is shown in Eig. [2](b), from which 
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FIG. 2. Edge spectrum and edge state SP angle 0 for mag¬ 
netically doped Sb 2 Te 3 QAH systems are plotted with: (1) 
y = 0 eV in (a) and (c). (2) E = 0.1 eV in (b) and (d). In 
(c) and (d), SP angle of edge mode I, II and III are plotted 
in red, blue and green. 


one can see that the number of left and right movers is 
unchanged, but their energy dispersions are shifted. 

After obtaining the edge state Green function, it is 
straightforward to calculate the corresponding SP vector 
S. First of all, Sx is vanishingly small compared to other 
spin components. This confirms the results from the ef¬ 
fective four band model. Thus, S only appears in the y-z 
plane and can be characterized by its magnitude |5'| and 
a polarization angle 0 relative to -\-z axis. For Sy, we no¬ 
tice that it vanishes at y = 0 for the effective four band 
model, while it is non-zero in the realistic tight-binding 
model. This difference is because of the out-of-plane mir¬ 
ror symmetry in our four band model, which is absent in 
realistic crystals. Therefore, both Sy and Sz are non¬ 
zero even at y = 0 in realistic tight-binding model. We 
plot Fermi energy Ep dependence of polarization angle 
0 at the local gate y = 0 eV and V = 0.1 eV in Fig. [ 2 ] 

(c) and (d). The polarization angles 0 for the modes I 
(red lines) and III (green lines) behave similarly to each 
other, while that of the mode II (blue lines) reveals a 
completely different characteristics. This is a clear ev¬ 
idence that the modes I and III are connected to each 
other, forming a non-chiral edge mode, while the mode 
II can be identified as the non-trivial GEM. In Fig. [2] 

(d) , we notice an abrupt change of polarization angle 0 
for the modes I and II when Ep between 0 meV and 
20 meV. Compared with the energy dispersion in Fig. [ 2 ] 
(b), we find that this change results from the strong hy¬ 
bridization effect between the modes I and II. Thus, the 
GEM and non-chiral mode are not well defined in this 
regime. In the other regime, we find a smooth change 
of SPs with respect to local gate voltages and chemical 


potentials. 

Experimental deteetion of edge spin polarization - Our 
theoretical calculations based on both the effective model 
and realistic tight-binding model have clearly demon¬ 
strate electrically tunable SP in magnetically doped 
(Bi,Sb) 2 Te 3 . However, the experimental detection of SP 
is not an easy task since the bulk is FM and it is un¬ 
clear how to distinguish SP of OEMs from that of bulk 
ferromagnetism by magnetization measurement. In con¬ 
trast, when the Fermi energy is tuned into the bulk gap, 
the transport signals are dominated by CEMs. Thus, it is 
more promising to search for SP of CEMs from transport 
measurements. 

Our proposal is based on a four terminal device with 
a FM probe as the fourth probe, as shown in Fig. [3] (a). 
This experimental setup is similar to that in the study 
of disordered leads in the QH system and here FM leads 
play the role of disordered leads [23|. When SP of the 
GEM is parallel to magnetization M of the FM lead, it 
can flow into the lead, while when they are opposite, it 
will be scattered. We apply a voltage drop between the 
leads Vi and Vs (Vi = V and ys = 0) and also introduce 
a split gate, denoted as SG in the Figl3] (a). Due to the 
split gate, there are two types of currents flowing from 
the lead y 2 to V 4 : the current which goes through V 3 , 
and the current ^2 flowing directly from^ to V 4 . Based 
on the Landauer-Biittiker formula 2J, [ 23 , the current 
ii to the lead V 4 shares the same chemical potential as 
ys = 0 , while the current ^2 is determined by the chem¬ 
ical potential in the lead V 2 = Vi = V. Importantly, 
we assume that chemical potentials of two currents H 
and ^2 have not reached equilibrium when they enter the 
FM lead. This is determined by the inelastic scatter¬ 
ing length, which is estimated as several jim for the QH 
case [ 261 , and we expect a similar length scale in our 
case. Thus, the corresponding SPs are also expected to 
be different for these two currents. Since the scattering 
rate into the FM lead V 4 depends on the relative angle 
between the SP of GEMs and M. We expect the trans¬ 
missions of y 2 ^ V 4 and Vs ^ V 4 also depend on M of 
the FM lead. As a result, the chemical potential in V 4 
will vary when rotating M. This provides a detectable 
signal, which is similar to the spin valve effect 0 , in 
transport measurements and can be directly related to 
SP of GEMs. It should be emphasized that the split gate 
SG plays an essential role here. Without the split gate, 
all currents flowing into V 4 come from Hs, and thus the 
corresponding chemical potential in V 4 must be equal to 
that of Vs , independent of magnetization direction of the 
FM lead. 

To formulate this idea, we assume that only p (p < 
1 ) fraction of the net edge current io can flow from y 2 
to Vs due to the split gate, as shown in Fig. [ 2 ] (a), so 
ii = pio^h = (1 — p)io. We further assume the q' (q) 
fraction of the current ii (^ 2 ) has spin polarizing along 
M. Since SP of GEMs depends on chemical potential, it 
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FIG. 3. Our proposed four terminal transport setup is shown 
in (a) with lines showing current flows. In (b), we show the 
spin valve effect at the FM lead V 4 : Only current with spin 
parallel to M can flow into V 4 , while those with spin anti¬ 
parallel to M will be reflected. 0 dependence of transmission 
qiS at Ef = 33 meV are plotted in (c). Based on this, we plot 
the evolution of transverse resistance R 24 with 0 in (d). 


is reasonable to assume q ^ q'. Therefore, the currents 
from V 2 (V 3 ) to V 4 are q'ii and gi 2 , respectively, as shown 
in Fig. [3] (a). The Hall conductance G 24 can be derived 
based on Landauer-Biittiker formula [l^ as 


G24 


h _ ^ (1 


(3) 


2 

with the conductance quanta Go = ^. From Eq. (|3)) , one 
can see that when there is no split gate (p = 1 ), G 24 = 

2 

Go = -I", which recovers the quantized Hall conductance 
and is independent of its spin information. When a split 
gate is introduced, 0 < p < 1 and q q'^ and G 24 will 
deviate from the quantized value and we discuss how to 
extract the information of SP of OEMs from the Hall 
resistance measurement in the Supplementary Materials 
0. 

Eor the realistic systems of magnetically doped 
(Bi,Sb) 2 Te 3 , we have shown additional non-chiral modes 
coexisting with OEMs. Thus, it is natural to ask how 
these non-chiral modes influence the above transport 
measurement. We consider one pair of non-chiral edge 
mode (mode I and HI in Eig. [2] (a)) and one GEM (mode 
H) and further assume no inter-channel scattering be¬ 
tween different modes for simplicity. We take the modes 
I and H flowing clockwise (as shown in Eig. [3] (a)), and 
the mode HI flowing counter-clockwise (flipping the di¬ 
rections of currents in Eig. [3] (a)). ^ 1 , 2,3 (^ 123 ) defined 
as the transmission into EM probe of edge current ^2 (g) 
for the modes I,H and HI, with Q = qi E q 2 ^ Q' = q'l E q 2 
and p = 1—p. The explicit expression of transverse resis¬ 
tance is shown in the Supplementary Materials [l3[. In 
the limit where the split gate vanishes (p = 1 — p = 1 ), 


the transverse conductance G 24 becomes 

7{Q'+ q3)-SqsQ' 

= ft- 20 ^, - 

which deviates from the quantized value. Thus, in con¬ 
trast to the single CME case. Hall resistance will depend 
on the magnetization direction of EM leads even without 
any split gate for the case with both CME and non-chiral 
modes. 

To apply Eq. (|1]) to magnetically doped Sb 2 Te 3 films, 
we need to extract the coefficients qiS from our realistic 
tight-binding model. As discussed above, g^s are deter¬ 
mined by the projection of SP of OEMs into the magneti¬ 
zation direction M of EM leads. Since SP of OEMs only 
exists in the y-z plane, we only concern the projection of 
SP into the y-z plane. Let’s assume M has an angle 0 
relative to + 2 ; direction in the y-z plane, as shown in Eig. 
[ 3 ] (b). The corresponding projection operator for OEMs 
is defined as: 


= I I, with I te) = I (5) 


Consequently, the angle dependent transmission qi (0) for 
the ith edge mode is given by 


qi{&) = 




(6) 


where lipi) is the wave function for the mode i. With 
Eq. m and Eq. [ 6 l we can calculate the transmissions 
qi{0)s for the modes I, H and HI as a function of 0 and 
Eermi energy for the local gate voltage V = 0.1 eE for 
our realistic tight-binding model [l3|. Eor the chemical 
potential Ep = 33 meV, 0 dependence of transmission 
g^s and the Hall resistance i? 24 , are shown in Eig. [3] (c) 
and (d), respectively. R 24 shows a strong dependence 
on thus revealing the spin valve effect for CEMs 0 . 
This provides a very clear and experimentally feasible 
evidence to detect spin signal in a QAH insulator. 

Discussions - In summary, we have shown that SP of 
CEMs can be generated, manipulated and detected in a 
QAH insulator in a fully electric manner. This paves the 
way of potential applications for the QAH effect in spin- 
tronics. Disorder is inevitable in realistic samples and 
we show the stability of SP of CEMs against disorder in 
the Supplementary Materials 0 . Although we focus on 
the magnetically doped (Bi,Sb) 2 Te 3 films here, the elec¬ 
tric controllability of SP of CEMs is a general property 
of a QAH insulator. The bulk topology (non-zero Chern 
number) only guarantees the existence of CEMs, but the 
detailed form of wave functions of CEMs depends on lo¬ 
cal electric environment, and thus can be controlled by 
a local gate voltage. We expect a similar effect also oc¬ 
curs in other QAH insulators, such as Mn doped HgTe 
quantum wells 0 and InAs/GaSb quantum wells 
where the local gate voltage can induce a local Rashba 
spin-orbit coupling and tilt SP of CEMs. 
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Analytical solution of edge wave function in effective model 

To solve for the zero mode wave function of the effective model, we can first apply a unitary transformation t/, 

10 1 0 \ 

0 10-1 

0 10 1 ’ 

10-1 0 / 



leading to 


H' = UHeffU^ 


rtik + Mz ivk- M- 0 \ 

—ivk-^ —rrik — Mz 0 M+ 

M+ 0 rrik — Mz —ivk-^ 

\ 0 M_ ivk- —ruk + Mz) 


We set M+ = M_ = 0 and H becomes block diagonal. We take m(k) = mo + m 2 k^ and assume mo ,2 > 0, and 
band inversion happens when \Mz\ > mo- We would like to emphasize that the two-band model taken in Ref. 
[l2| is just one block of the above four-band Hamiltonian ([7|). Nevertheless, the bases for the Hamiltonian © are 
(|+, t), | —,1), |-h,i), | —, t)): where ” denotes bonding/anti-bonding state between top and bottom surface states. 
The corresponding spin operators are given by 



/O 0 1 0\ 
0 0 0 1 


/O 0 -i 0\ 
0 0 0 z 


/l 0 0 0\ 

0-100 

— 

10 0 0 
1^0 1 0 Oyi 

5 Sy = 

i 0 0 0 

V -i 0 oj 

= 

0 0-10 
VO 0 1/ 


Thus , we find that the spin operators in the four band model are different from those used in the two-band model in 
Ref. [l2|, which are more relevant to the cold atom system. 

Now we consider the case with Mz > 0, in which the lower block of H' is in the QAH phase while the upper block 
is in the trivial phase. For the lower block, we have 


fruk- Mz -ivk^ A 
ivk- -rrik + Mz) 


{rrik - Mz)(Jz + v{kyax -f k^cTy). 


We choose the following set of parameters which satisfies the above band inversion condition: u = 1, Mz = 0.5, mo = 
0.1, m 2 = 0.25. 

The wave-function of zero energy mode of the chiral edge mode can be solved for the lower block Hamiltonian. Let 
us assume the edge is along the y direction, so we consider kx = 0 first, and check the eigen equation for zero mode 

[(mo - m 2 dl - Mz)(Tz - ivxCTxdyjtpiy) = 0 (7) 

Multiply both sides by 

[(mo - m 2 dy - Mz)ay + Vz;dy\xl){y) = 0 (8) 

Consider the ansatz ^lr{y) = where ay^s = s = ±1, we have 

A± = ;;^(-sv ± x/v2 +4m2(mo - Mz)) (9) 

2m2 

Thus the wave function should be the following linear combination: 

^(y) = '^ica+e~^+^ + Cs-e~^-y)^s 


( 10 ) 
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©(degree) ©(degree) 

FIG. 4. In (a), we plot 0 dependence of transmission q at various chemical potential for the chiral edge state in the effective 
model. In (b), we fix the chemical potential and plot the 0 dependence of q. Transverse resistance R 24 at different split gate 
transmission p is plotted in (c). 


with 


A^ + Ai = 
A^Ai = 


sv 

m 2 

mo - Mz 
m 2 


( 11 ) 


Assuming the material (vacuum) region is at x > 0 (x < 0), then we require that wave function vanishes at both 
X = +00 and X = 0. This requires Cs+ = —Cs_, Aj- > 0 and s to be negative. So the wave function is given by 


= c(e“^+^ — (12) 

Notice that the spinor part is defined in the hybridized bases (| + i), | — t))- Then can be written as 

'^1+ i) + 1“ t) = '^(1^ i) + 1^ i)) + (1^ t) ~ 1^ t)) 

= \t) (g) {i\ i) + I t)) - 1^) G) (-^1 I) + I t)) 

= i|t)(8)(| ty))+i|^)G)(|4)) (13) 

We emphasize again that the above solution does not correspond to a net spin polarization along the y direction, as 
discussed in Ref. Hi , due to the different bases. The wave function of OEMs consists of two parts: one part on the 
top surface with SP along the -^y direction and the other on the bottom surface with SP along the —y direction. 
These two contributions exactly cancel each other, leading to zero net SP of OEMs. Erom the symmetry perspective, 
we notice that there is an emergent mirror symmetry between the top and bottom surfaces, which requires Sy = ^. 
Therefore, applying a local gate voltage breaks this mirror symmetry and yields an imbalance of SP between the top 
and bottom surfaces, giving rise to a finite Sy. The absence of Sz is completely accidental, as a result of setting 
Ao,i = h = 0. Numerically, when assigning non-zero values to Ao,i and /i, non-zero Sz appears, as shown in Eig. [ 1 ] 
in the main text. Therefore, Sz is insensitive to this mirror symmetry breaking, just as we have numerically verified. 
To conclude, applying a local gate voltage changes both magnitude and direction of GEM SP. 


TRANSVERSE CONDUCTANCE FROM LANDAUER-BiiTTIKER FORMULA 


In this section, we will derive analytical expressions of transverse conductance based on Landauer-Buttiker formula 
3,0. We will first discuss the situation with one single chiral edge mode, and then generalize the expression to a 
case where a chiral edge mode and one pair of non-chiral edge mode coexist. 

Following the setup configuration in Fig. [3] (a) in the main text, the conductance matrix G for a single chiral edge 
state (flowing clockwise) is given by 


/O {l-p){l-q) p{l-q') il-p)q + pq'\ 

10 0 0 

Op 0 0 

\0 (1 - p)q pq' 0 / 


G = Gqx 
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where the bases for matrix G is (Vi, V 25 V 3 , ^ 4 ). Here p is the split gate transmission, while q' and q are transmissions 
into FM probe V 4 of currents from V 2 and V 3 
from probe Vj to probe Vi. Then, current voltage relations are 


Here a matrix element Gij describes the conductance contribution 


AA 

I 2 

h 

\h) 


(I -p{l-q') -{l-p)q-pq'\ 


— Go X 


-1 

0 

VO 


1 

-p 

-{l-p)q 


0 

P 

-pq' 


{l-p)q+pq' / 


(Vi\ 

V 2 

\vp 


Here U {i G 1, 2, 3,4) denotes the current that flows into the ith probe. Notice that /2 = /4 = V 3 = 0 and h = 
giving rise to 

^1 = ^2 

(1 -p)qV 2 = [(1 -p)q+pq']V4 
h = GopVi 


-h, 


The transverse conductance 6 ^ 24 ”^*^* is given by 


G chiral _ 

24 — 


h 


Vo-Va 


= Go 


{l-p)q+pq' 

q' 


(14) 


(15) 


Numerically, we plot the 0 dependence of transverse resistance i ?24 in Fig. |4] (c). Here, we take current ii (from 
V 3 ) at Ep = 0.1 eV and 12 (from V 2 ) at Ep = 0.4 eV. Transmission of ii and 12 into V 4 is calculated in Fig. |4] (a), 
which will be discussed in details in a later section. In the Fig. |4] (b), we extract the transimission information for 
Ep = O.leV (blue line) and Ep = 0.4eV (red line) from Fig. |l](a), which provides us with the q and q' at different 0. 
Transmission into the split gate p is chosen at various values. As long as the split gate structure exists (p 7 ^ 0), there 
will be an obvious 0 dependence of Hall resistance i? 24 , as shown in Fig. [4| (c). But when the split gate structure is 
absent, transverse resistance will be quantized despite the value of 0. 

For the magnetically doped Sb 2 Te 3 QAH system, we consider two clock-wise modes (mode I and HI) and one 
counter-clockwise mode (mode II). Let us define gi, 2,3 (^ 1 , 2 , 3 ) transmission into FM probe of edge current 

12 (H) for the modes I,H and HI, with Q = qi E q 2 ^ Q' = Qi + q2 and p = 1 — p. So the conductance matrix for mode 
I and mode II can simply take the form of Eq[T4j while conductance matrix for mode HI is as follows: 


G/ = Go X 


/O (l-p)(l-gi) pil-q'i) (1 
1 0 0 
Op 0 

\0 (1 - p)qi pq{ 

Gii = Gi{qi q2, q'l ^ q' 2 ) 

( 0 1 0 0 \ 
(1 -P)(l - 93 ) 0 p {l-p)qz 
p{l -qs) 0 0 pq 3 
qs 0 0 0 


■p)qi +pq'i\ 

0 


J 


Gill = Go X 




J 


The net conductance matrix is simply an addition of all three conductance matrix Gnet = Gj E Gn E Gm. After 
some calculations, we arrive at the following Hall conductance 


g; 


non-chiral 

24 


(7 + 2 p) [Q' -j- p{Q — Q')] + ^3 [7 - 1 - 2 p — 3Q' E p( 2 p E I) (Q^ “ Q)] 


h (2E p)Q' — ^ 3(1 + pQ') 

In the limit where the split gate vanishes (p = I — p = I), the transverse conductance (^non-chirai becomes 


G non-chiral 
9 


7((3' + 93 ) - iqoQ' 


^24 


h 


2Q' - qo 


(16) 


(17) 


ELECTRON CURRENT TRANSMISSION AT THE FERROMAGNETIC PROBE 


In the previous discussion, we have shown that transverse conductance G 24 has a strong dependence on transmissions 
at the FM probe V 4 . In the main text, we have defined the transmission qi based on a projection operator approach. 
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Mode2 
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0(degree) 


FIG. 5. Chemical potential and 0 dependence of transmissions of mode I, II and III are plotted in (a), (b) and (c). Here we 
have applied a local gate V = 0.1 eV throughout the calculations. In (d), we compare the transmission of all three modes at a 
fixed Fermi energy Ef = S3 meV. 


Since the edge state Green function can be calculated with the help of iterative Green function method, we can easily 
rewrite the transmission expression in terms of the edge state Green function Gi{uj) for edge mode i (i=I,II,III) 


qi{S) 




(18) 


Based on Eq. [181 Fig. |l](a), we calculate the chemical potential and 0 dependence of transmission for the chiral 
edge state in the effective QAH model. Here we apply a local gate voltage V = 0.3 eV. The blue line {Ep = 0.1 
eV) and red line {Ep = 0.3 eV) denote two different chemical potentials due to the split gate effect. The two Fermi 
energies are the Eps we choose to calculate transverse resistance in Fig|l](c). We further plot the transmissions q{Ep) 
at these two fixed Fermi energies in Fig. |4](a). 

For the magnetically doped Sb 2 Te 3 QAH system, non-chiral edge modes show up and we consider the co-existence 
of a pair of non-chiral edge modes (mode I and mode III) and a chiral mode (mode II). Following a similar approach, 
we make use of the edge state Green function and plot chemical potential and 0 dependence of transmissions g^s for all 
three modes in Figl5] (a) to (c). Here the white dashed line shows the fixed chemical potential Ep = 33 meV, which 
we picked to calculate i ?24 in the main text. In Figl5](d), we compare the transmission of all three modes at the fixed 
Fermi energy Ep = 33 meV. We notice that transmission of mode II peaks at around ^ = 100 degree while that of 
mode I and III peaks at around 0 = 200 degree. This means that transmission peaks happen at the polarization angle 
of these edge modes. The above results are consistent with the fact that FM probe is transparent to the electrons 
with spin parallel to the magnetization of the FM probe. 


MAGNITUDE OF EDGE SPIN POLARIZATION FOR MAGNETICALLY DOPED SB2TE3 QAH 

SYSTEMS 

We have found that for magnetically doped Sb 2 Te 3 QAH systems, spin polarization S is a two dimensional vector 
living in y-z plane, due to the absence of its x component. Therefore, S can be written in terms of polar coordinates 
S = (|5'| cos6>, |5'| sin6>). In the main text, we have plotted the polarization angle 0 of mode I, II and III. It is shown 
that mode I (left mover) and mode III (right mover) share a very similar 0 characteristics, and we claim that these 
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Ep(meV) E/r(meV) 

FIG. 6. Amplitude \S\ of spin polarization for magnetically doped Sb 2 Te 3 QAH systems with V = 0 eV and y = 0.1 eV are 
shown in (a) and (b). \S\ of edge mode I, II and III are plotted in red, blue and green. 



FIG. 7. In (a), the density of states of the nth lattice sites pn is plotted in units of pi, for n = 1, 2, 3,4, 5. It is easy to see that 
the decay length is simply one lattice constant. In (b), we compare the spin angle 0 of the n = 1 lattice sites (red line) and 
that of the entire contribution of n = 1, 2, 3,4, 5 lattice sites (blue line). 


two modes should be connected to each other and become trivialized. To make our results complete, here we show 
that magnitude plot of spin polarization in Fig. [6l Following the same color convention, \S\ of edge mode I, II and 
III are plotted in red, blue and green. In fact, the value of \S\ can be closely related to the calculation details. For 
example, when locating the OEMs in the iterative Green function (DOS) calculation, a small error in locating the 
peak position of DOS may lead to a different \S\ value of the corresponding GEM. Compared to jEj, polarization 
angle 0 is barely affected by such calculation error. Therefore, we would like to point out that polarization angle 0 is 
more essential in our discussion here. 


DECAY LENGTH OF CHIRAL EDGE STATES 

Let us denote the density of state calculated from the nth lattice sites as with n = l,2,3,.... In FigI3 (a), we 
have plotted the value of pn in units of pi of mode II in the Sb 2 Te 3 based QAHE system. It clearly shows that the 
edge state is highly localized on the edge, and exponentially decays into the bulk. By doing an exponential fit, we find 
that the decay length is 1.03 times lattice constant. In our calculation for SP, we have included the total contribution 
from all n = 1,2, 3,4, 5 lattice sites, which we identify as the major contribution. In general one can integrate over the 
entire semi-infinite space, which is rather costly in calculation, while giving only a small modification to our present 
results. Actually, we can even show the SP angle 0 is dominated by the contribution from the n = 1 lattice sites. As 
shown in Fig. [7| (b), we calculated the SP angle 0 of CEM (mode II) with Vgate = 0.1 eV for only n = 1, and plot 0 
evolving with Fermi energy (red line). In comparison, we also plot the SP angle calculated for all n = 1, 2, 3,4, 5 lattice 
sites with blue line. It is clear that the n = 1 lattice sites give the major contributions to 6>, while n = 2,3,4, 5 lattice 
sites together only give a rather small correction to 0. We also would like to emphasize that we are not interested the 
exact value of CEM SP, which depends on the details of samples in experiments. The essential physics here is that 
the CEM SP can be manipulated electrically, despite its exact value. 
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FIG. 8. We plotted possible reconstruction of chiral edge states under the influence of FM probe. The yellow lines are the 
chiral edge states. When FM probe destroy the QAHE locally, a new boundary emerges and edge states appear on this new 
boundary. 


INFLUENCE OF FM PROBE ON THE QAHE 

It is possible that attaching a FM gate to a QAH sample may change the QAH sample locally. For example, the 
value of CEM spin polarization may be affected in this regime. However, what is physically interesting is the tunability 
of spin polarization as a result of electrical approaches (chemical potential and local gate tuning). The exact value of 
CEM spin polarization is not that interesting to us. Therefore, we expect the angle-dependent transverse resistance 
should still be observed in this case. 

In an extreme situation, QAH state is destroyed locally in the vicinity of FM probe, as shown in Fig. [HI In this case, 
one can simply treat this local part of QAH sample as an extension of the FM gate. Then a new boundary between 
FM probe (including its extension) and the QAH sample emerges and edge states still exist on this new boundary, as 
long as the rest of sample remains in the QAH phase. As a result, our physical picture of this transport experiment 
still holds. 


STABILITY OF EDGE SPIN POLARIZATION UNDER DISORDER 

In this section, we show that the above spin polarization of QAH edge states in the simple effective model (Eq. 
[1] in the main text) is robust under disorder effect. We start from a finite size real space Green function with a 
N, X Ny lattice configuration with = 50 and Ny = 2. To discuss its chiral edge states, we perform iteration 
of Green function along the y direction, so that system is finite along the x direction while semi-infinite along y 
direction. In this case, SDOS (spin polarization) of edge states on the ^-edge will be obtained. Disorder potentials 
are white-noise-like on-site potentials, and are included locally on the edge for simplicity. The magnitude of disorder 
potential is between —W and VF, where W is the disorder strength. We also perform the configuration averaging, the 
number of which is defined as Nave, foi* disorders. 

In Fig. [9] (a), we first calculate edge spin polarization without any disorder {W = 0). A finite local gate voltage 
{Vgate = 0.1 cV) and a finite Fermi energy {Ep = 0.1 eV) are considered to guarantee non-zero spin polarization. The 
x-axis represents lattice sites of the one-dimensional lattice chain on the ^-edge and ^-axis is for spin polarization. 
We can clearly see strong finite size effect on the ID edge chain: At both ends of this ID chain, there are abrupt 
changes of spin polarization in all three components. Such a finite size effect can be removed by imposing periodic 
boundary condition along x direction, and resulting spin polarization is exactly shown in Fig. [1] (a) and (b) (in the 
main text). This means that for a finite size QAH sample, spin polarization of an electron on the ^-edge is along -^y 
direction. As shown in Fig. [TOl while approaching the sample corner between ^-edge and x-edge, its spin polarization 
(green arrows) gradually becomes along — z direction. Finally, when this electron enters x-edge, its spin polarization 
becomes along —x direction. 

When disorder is introduced, we can see that the fluctuation of spin polarization increases as the disorder strength W 
increases, as shown in Fig. |9](b)-(e). But this disorder induced fluctuation decreases as we increases the configuration 
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FIG. 9. Spin polarization of effective QAH model with various disorder strengths W and Nave = 100 are plotted: (a) VF = 0 
eV, (b) VF = 0.1 eV, (c) IF = 0.2 eV, (d) IF = 0.3 eV, (e) IF = 0.4 eV. In (f), we increases Nave to 1000 with disorder strength 
IF = 0.4 eV. 



FIG. 10. In this figure, we give a systematic plot of the edge spin polarization of a finite size QAH sample based on Fig. [9] (a). 
Red arrow represents edge currents that flow around the QAH sample, while green arrows represent the corresponding spin 
polarization vectors of the edge current. On the right hand side of the QAH sample, we show in details how the SP of a chiral 
edge current on the y edge evolves due to the finite size effect. 


averaging number Nave (Fig. [9] (f)). This concludes that spin polarization of chiral edge state in QAHE is robust 
against disorder effects. 
























